Maximal entanglement versus entropy for mixed quantum states 
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Maximally entangled mixed states are those states that, for a given mixedness, achieve the greatest 
possible entanglement. For two-qubit systems and for various combinations of entanglement and 
mixedness measures, the form of the corresponding maximally entangled mixed states is determined 
primarily analytically. As measures of entanglement, we consider entanglement of formation, relative 
entropy of entanglement, and negativity; as measures of mixedness, we consider linear and von 
Neumann entropies. We show that the forms of the maximally entangled mixed states can vary 
with the combination of (entanglement and mixedness) measures chosen. Moreover, for certain 
combinations, the forms of the maximally entangled mixed states can change discontinuously at a 
specific value of the entropy. 
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I. INTRODUCTION 

Over the past decade, the physical characteristic of the 
entanglement of quantum-mechanical states, both pure 
and mixed, has been recognized as a central resource 
in various aspects of quantum information processing. 
Significant settings include quantum communication jl), 
cryptography teleportation Q, and, to an extent that 
is not quite so clear, quantum computation [Q. Given the 
central status of entanglement, the task of quantifying 
the degree to which a state is entangled is important for 
quantum information processing and, correspondingly, 
several measures of it have been proposed. These in- 
clude: entanglement of formation |3|, |6[, entanglement of 
distillatio n \f\ , relative entropy of entanglement ^ , neg- 
ativity H [uJL and so on. It is worth remarking that even 
for the smallest Hilbert space capable of exhibiting entan- 
glement, i.e., the two-qubit system (for which Wootters 
has determined the entanglement of formation ||), there 
are aspects of entanglement that remain to be explored. 

Among the family of mixed quantum mechanical 
states, special status should be accorded to those that, 
for a given value of the entrop y |fl3T| , have the largest pos- 
sible degree of entanglement |12]j. The reason for this is 
that such states can be regarded as mixed-state gener- 
alizations of the Bell states, the latter being known to 
be the maximally entangled 2-qubit pure states. The 
notion of maximally entangled mixed states was intro- 
duced by Ishizaka and Hiroshima jl3| in a closely related 
setting, i.e., that of 2-qubit mixed states whose entan- 
glement is maximized at fixed eigenvalues of the density 
matrix (rather than at fixed entropy of the density ma- 
trix). Evidently, the entanglement of the maximally en- 
tangled mixed states of Ishizaka and Hiroshima cannot be 
increased by any global unitary transformation. For these 
states, it was shown by Verstraete et al. [|l4| that the max- 
imally property continues to hold if any of the following 



three measures of entanglement — entanglement of forma- 
tion, negativity, and relative entropy of entanglement — is 
replaced by one of the other two. 

The question of the ordering of entanglement measures 
was raised by Eisert and Plenio Jl^] and investigated 
numerically by them and analytically by Verstraete et 
al. iQ. It was proved by Virmani and Plenio |l7| that 
all good asymptotic entanglement measures are either 
identical or fail to uniformly give consistent orderings of 
density matrices. This implies that the resulting maxi- 
mally entangled mixed states (MEMS) may depend on 
the measures one uses to quantify entanglement. More- 
over, in finding the form of MEMS, one needs to quantify 
the mixedness of a state, and there can also be ordering 
problems for mixedness. This implies that the MEMS 
may depend on the measures of mixedness as well. 

This Paper is organized as follows. We begin, in 
Sees. II and III, by reviewing several measures of en- 
tanglement and mixedness. In the main part of the Pa- 
per, Sec. 



IV 



we consider various entanglement-versus- 
mixedness planes, in which entanglement and mixedness 
are quantified in several ways. Our primary objective, 
then, is to determine the frontiers, i.e., the boundaries 
of the regions occupied by physically allowed states in 
these planes, and to identify the structure of these max- 
imally entangled mixed states. In Sec. [v| we make some 
concluding remarks. 



II. ENTANGLEMENT CRITERIA AND THEIR 
MEASURES 

It is well known that there are a large number of 
entanglement measures E. For a state described by the 
density matrix p a good entanglement measure must 
satisfy, at least, the following conditions jl8], fll| : 
CI. (a) E{p) > 0; (b) E(p) = if p is not entangled M; 
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and (c) E^Bell states) = 1. 

C2. For any state p and any local unitary transformation, 
i.e., a unitary transformation of the form U a ® Ub, the 
entanglement remains unchanged. 

C3. Local operations, classical communication and 
postselection cannot increase the expectation value of 
the entanglement. 

C4. Entanglement is convex under discarding informa- 
tion: J2iPi E (pi) ^ E (J2iPiPi)- 

The entanglement quantities chosen by us satisfy the 
properties C1-C4. Here, we do not impose the condition 
that any good entanglement measure should reduce to 
the entropy of entanglement (to be defined in the follow- 
ing) for pure states. 



A. Entanglement of formation and entanglement 
cost 

The first measure we shall consider is the entanglement 
of formation Ep B ; it quantifies the amount of entangle- 
ment necessary to create the entangled state. It is defined 
by 



ip)= min V ftj B(|^){^|), (2.1) 



where the minimization is taken over those probabilities 
{Pi} and pure states {ipi} that, taken together, reproduce 
the density matrix p — YliPi\' , Pi)(''Pi\- Furthermore, the 
quantity E(\tpi)(ipi\) (usually called the entropy of en- 
tanglement) measures the entanglement of the pure state 
\ipi) and is defined to be the von Neumann entropy of the 

reduced density matrix p[ A ^ = TxelV'iKV'il) i-e., 

S(|^)(^l) = -Trp| A) log 2 ^. (2.2) 

For two-qubit systems, Ep can be expressed explicitly 
as § 

MP) = h Q[l + y/1 - C(p)*]) , (2.3a) 
h(x) = — xlog 2 x — (1 — x) log 2 (l — x), (2.3b) 
where C(p), the concurrence of the state /?, is defined as 
C(p) = max{0, \f\~i- \/% - \/% - a/a7}, ( 2 -3c) 



in which Ai,... , A4 are the eigenvalues of the matrix 
p(o-y 8> a y )p*(o~ y ® o~y) in nondecreasing order and a y 
is a Pauli spin matrix. Ep(p), C(p), and the tangle 
r(p) = C(p) 2 are equivalent measures of entanglement, 
inasmuch as they are monotonic functions of one another. 

A measure associated with the entanglement of forma- 
tion is the entanglement cost Eq [^) which is defined via 



E c {p) = lim 



Ep{p® n ) 



(2.4) 



B. Entanglement of distillation and relative 
entropy of entanglement 



Related to the entanglement of formation is the en- 
tanglement of distillation Ep, H , which characterizes the 
amount of entanglement of a state p as the fraction of 
Bell states that can be distilled using the optimal purifi- 
cation procedure: Ep,(p) = linin^oo m/n, where n is the 
number of copies of p used and m is the maximal num- 
ber of Bell states that can be distilled from them. The 
difference E-p — Ep, can be regarded as undistillable en- 
tanglement. Ep, is a difficult quantity to calculate, but 
the relative entropy of entanglement £"r |^|, which we 
shall define shortly, provides an upper bound on Ep> and 
is more readily calculable than it. For this reason, it is 
the second measure that we consider in this Paper. It is 
defined variationally via 



E R (p) = minTr (plogp - ploger) , (2.5) 



where D represents the (convex) set of all separable den- 
sity operators a. In certain ways, the relative entropy of 
entanglement can be viewed as a distance T>(p\\o~*) from 
the entangled state p to the closest separable state a*. 
We remark that for pure states of two-qubit systems the 
relative entropy has the same value as the entanglement 
of formation. 



C. Negativity 

The third measure that we shall consider is the negativ- 
ity. The concept of the negativity of a state is closely re- 
lated to the well-known Peres-Horodecki condition for the 
separability of a state pl[ . If a state is separable (i.e., not 
entangled) then the partial transpose of its density ma- 
trix is again a valid state, i.e., it is positive semi-definite. 
It turns out that the partial transpose of a non-separable 
state has one or more negative eigenvalues. The nega- 
tivity of a state indicates the extent to which a state 
violates the positive partial transpose separability crite- 
rion. We will adopt the definition of negativity as twice 
the absolute value of the sum of the negative eigenvalues: 



N(p) = 2max(0,-A neg ), 



(2.6) 



This is the asymptotic value of the average entanglement 
of formation. Eq is, in general, difficult to calculate. 



where A nog is the sum of the negative eigenvalues of p TB . 
In C 2 ® C 2 (i.e., two-qbit) systems it can be shown that 
the partial transpose of the density matrix can have at 
most one negative eigenvalue (see App. |X|). It was proved 
by Vidal and Werner |l(| that negativity is an entan- 
glement monotone, i.e., it satisfies criteria C1-C4 and, 
hence, is a good entanglement measure. We remark that 
for two-qubit pure states the negativity gives the same 
value as the concurrence does. 
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D. Bures metric 

The Bures metric of entanglement is defined as 

E B (p) = min(2-2 y /F(p,a)), (2.7) 

where F{p,a) = (Tr \J ^fap^fa) 2 is the fidelity. In the 
same way that relative entropy can, this entanglement 
measure can be viewed as the distance from the clos- 
est separable state to the entangled state considered, 
where the distance is now defined by 2?(p||er) = (2 — 
2y/F(p, ct)) [ pi) . We remark that for two-qubit pure 
states the Bures metric reduces to the tangle defined in 
Sec. HA. 



E. Lewenstein-Sanpera entanglement 




FIG. 1: Comparison of negativity and concurrence. Dots rep- 
resent randomly generated states. (The form of the boundary 
is discussed in the text.) It is apparent that these two mea- 
sures of entanglement can give different orderings for pairs of 
states. 



It was shown by Lewenstein and Sanpera fl23fl that any 
density matrix p has a decomposition into two parts: 



p = \p s + (I - X)p c , 



(2.8) 



where p s is separable, p c is entangled, and the weight A 
is maximal, in which case the decomposition is unique. 
They refer to p s as the best separable approximation 
(BSA) to p. It should be pointed out that, in general, it 
is nontrivial to establish the decomposition, even in the 
simplest relevant setting of C 2 ® C 2 systems. Evidently, 
A and p c contain information about the entanglement of 
p. Karnas and Lewenstein later showed that the 
quantity E'ls = (1 — A), which we will call the LS entan- 
glement, satifies the above criteria, and hence is a good 
entanglement measure. For C 2 ® C 2 systems it turns out 
that p e is a pure state, i.e., | , 0e)(V'e|, and in this case it 
suggests that the quantity £ls, defined via 



LS 



(l-A^dVOWel), 



(2.9) 



may also be a good entanglement measure. We remark 
that the for two-qubit case the entanglement measure 
(1 — A) is known to be equal to the Schmidt measure 
introduced in Ref. pE| . 

Even though the LS decomposition is not , in general, 
straightforward to find for the states in Eq. (L5), the LS 
decomposition reads 




fx + VabOO Vab \ 
aO 
0_ b 
ab y+VabJ 



+ (1-X)p e 



(2.10) 



where p c = \<f>+){<f>+\ (with |«£±) = (|00> ± |11»/V2), 
(1 — A) = r — 2\fab, and we consider only r — 2\fab > (as 
for r — 2\fab < the whole density matrix is separ able) . 

If we compute the concurrence of the state in Eq. (4.5), 
we find C — max{0,r — 2y/ab}. It is interesting to note 



that the LS entanglement gives, for the Ansatz states in 
Eq. (4.5), the same value as the concurrence does; but 



for an arbitrary state this is not, in general, true |26 



F. Ordering difficulties with entanglement 
measures 

We now pause to touch on certain difficulties posed 
by the task of ordering physical states using entangle- 
ment. As first discussed and explored numerically by 
Eisert and Plcnio 1 15 , and subsequently investigated an- 
alytically by Verstraete et al. |[(|, different entanglement 
measures can give different orderings for pairs of mixed 
states. This can be seen, e.g\, from the plot of concur- 
rence versus negativity, Fig. |l|. The upper boundary is 
readily seen to be N < C whereas the lower boundary 

can be derived, giving N > \J^{C~\) 2 + \ + (C-l); see 
Ref. |l(| . Hence, when we wish to explain maximally en- 
tangled mixed states we need to be very explicit about 
the measure of entanglement (and also mixedness; see the 
following section). Different measures are likely to lead 
to different classes of MEMS states. 

We end this section by mentioning the three entan- 
glement measures that we shall use to compute the 
entanglement-mixcdncss frontiers: entanglement of for- 
mation, negativity, and relative entropy of entanglement. 
The first two of these are straightforward to compute, at 
least in two-qubit settings. For the third, certain re- 
sults are available jl4|, |lj] that ease the computation of 
MEMS. We have also reviewed four additional measures 
(entanglement cost, entanglement of distillation, LS en- 
tanglement, and the Bures metric). Of these, however, 
the first two are rather difficult to compute, let alone 
maximize; the third is also difficult to compute, at least 
in practice. As for the fourth, calculating the entangle- 
ment involves finding the closest separable states (as is 
required for the case of relative entropy). 
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III. MEASURES OF MIXEDNESS 

In the entanglement-measure literature, two measures 
of mixedness have basically been used: (l — Tr[p 2 ]) and 
the von Neumann entropy. Whereas the latter has an 
natural significance stemming from its connections with 
statistical physics and information theory, the former is 
substantially easier to calculate. Of course, for density 
matrices that are almost completely mixed, the two mea- 
sures show the same trend. 



A. von Neumann entropy 

The von Neumann entropy, the standard measure of 
randomness of a statistical ensemble described by a den- 
sity matrix, is defined by 

Sv(p) = -Tr(plogp) = -J2>H log Xi , (3.1) 
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FIG. 2: Comparison of linear entropy and von Neumann en- 
tropy. Dots represent randomly generated states; pure (Rank- 
1) states lie at the origin; Rank-2 states lie on segment b; the 
lighter dots in the interior are Rank-3 states; the darker ones 
are Rank-4 states. The lower boundary comprises three seg- 
ments meeting at cusps, whereas the upper boundary is a 
smooth curve. The two dashed lines represent thresholds of 
entropies beyond which no states contain entanglement. 



where Xi are the eigenvalues of the density matrix p and 
the log is taken to base Af, the dimension of the Hilbcrt 
space in question. It is straightforward to show that the 
extremal values of Sy are zero (for pure states) and unity 
(for completely mixed states) . To compute the von Neu- 
mann entropy it is necessary to have the full knowledge 
of the eigenvalue spectrum. 

As we shall mention in the following subsection, there 
is a linear entropy threshold above which all states are 
separable. Qualitatively identical behavior is encoun- 
tered for the von Neum ann entropy. In particular, as we 
shall see in Sec. [V C 1 , for two-qubit systems all states 
are separable for 5 V > -(1/2) log 4 (l/12) w 0.896. 



B. Purity and linear entropy 

The second measure that we shall consider is called 
the linear entropy and is based on the purity of a state, 
V = Tr [p 2 ] , which ranges from 1 (for a pure state) to 
1/Af for a completely mixed state with dimension Af. 
The linear entropy Sh is defined via 



Af 



Af-1 



U-Trb 2 ]), 



(3.2) 



which ranges from (for a pure state) to 1 (for a max- 
imally mixed state). The linear entropy is generally a 
simpler quantity to calculate than the von Neumann en- 
tropy as there is no need for diagonalization. For C 2 (g> C 2 
systems the linear entropy can be written explicity as 



S L (p)^-(l-Tr[p 2 ]). 



(3.3) 



A related measure, which we shall not use in this Pa- 
per (but mention for the sake of completeness), is the 
inverse participation ratio. Defined via R = 1/Tr [p 2 ] , 



it ranges from 1 (for a pure state) to Af (for the maxi- 
mally mixed state). An attractive property of the inverse 
participation ratio is that all states with R > J\f — 1 arc 
separable ||, which implies all states with a linear en- 
tropy S L (p) > Af(Af — 2)/ (jV — l) 2 (which is 8/9 when 
Af = 4) are separable. 



C. Comparing linear and von Neumann entropies 

The aim of this subsection is to illustrate the differ- 
ence between the linear and von Neumann entropies. We 
shall do this by considering the Af = 4 Hilbert space, 
and seeking the highest and lowest von Neumann en- 
tropies consistent with a given value of linear entropy. 
Before restricting Af to 4, the corresponding stationarity 
problem reads: 



5[Sv(j>)+P 



AT -I 
2Af 



S h (p) -(v- l)Trp = 0, (3.4) 



where f3 and v are, respectively, Langrange multipliers 
that enforce the constraints that linear entropy be fixed 
and that p be normalized. Thus, we arrive at the engag- 
ing self-consistency condition 



p = cxp(-^ - f3p), 



(3.5) 



in which v and (1 can be fixed upon implementing the 
constraints. 

By hypothesizing certain relationships among the 
eigenvalues of the density matrix, we have been able, 
for the case of Af = 4, to find certain solutions of the 
self-consistency condition and, hence, a candidate for the 
boundary of the region of Sh vs. Sy plane (see Fig. ||) 
corresponding to legitimate density matrices. These so- 
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lutions, when given in terms of eigenvalues, read 



r 1 — r 1 — r 1—r 1 + 3r 

{r,l - r,0,0}, 
1—r 1 — r 



j, for 0< r < 1, (3.6a) 



for - < r < 1, (3.6b) 



for < r < -, (3.6c) 



f 1—r 1 — r ~i 
{^— — °}< 
( A—r 4— r 4— r 3r~i 

{■W^^}' fOT < r < 1, (3.6d) 



and they correspond to the upper boundary, and the low- 
est, middle, and highest pieces of the lower boundary, 
respectively. Note that the lower boundary comprises 
three (in general, AT — 1) segments that meet at cusps. 
We remark, parenthetically, that the solutions with zero 
eigenvalues correspond to extrema within some subspace 
spanned by those eigenvectors with nonzero eigenvalues, 
and therefore only obey the stationarity condition ( |3.5| ) 
within the subspace. 

Is there any significance to the boundary states? 
Boundary segment (a) includes the Werner states de- 
fined in Eq. fl4.7p . Boundary segment (b) includes the 
first branch of the MEMS for Ep and Sl specified below 
in Eq. (p~6|). The segment (c) includes the states 

Pc = r\<f> + ){<t>+\ + i^(|01)(01| + |10)(10|). (3.7) 

States on segment (d) are all unentangled. Of course, the 
boundary segments include not only the specified states 
but also all states derivable from them by global unitary 
transformation. 

As for the interior, we have obtained this numerically 
by constructing a large number of random sets of eigen- 
values of legitimate density matrices, and computing for 
each the two entropies. As Fig. [2] shows, no points lie 
outside the boundary curve, providing confirmatory evi- 
dence for the forms given in Eq. ( |3 ,6[ ) . 

The fact that the bounded region is two-dimensional 
indicates the lack of precision with which the linear en- 
tropy characterizes the von Neumann entropy (and vice 
versa, if one wishes). In particular, the figure reveals an 
ordering difficulty: pairs of states, A and B, exist for 
which S£ ~ and ~ differ in si § n - Worse still, 
states having a common value of Sy have a continuum 
of values of 5l, and vice versa. 



IV. ENTANGLEMENT- VERSUS-MIXEDNESS 
FRONTIERS 

We now attempt to identify regions in the plane 
spanned by entanglement and mixedness that are inhab- 
ited by physical states (i.e., characterized by legitimate 
density matrices) . We shall consider the various measures 
of entanglement and mixedness discussed in the previous 
section. Of particular interest will be the structure of 
the states that inhabit the frontier, i.e., the boundary 



delimiting the region of physical states. Frontier states 
are maximal in the following sense: for a given value 
of mixedness they are maximally entangled; for a given 
value of entanglement they are maximally mixed. 



A. Parametrization of maximal states 

The aim of this subsection is to deriv e the general form 
of the maximal states given in Eq. (4.4), which is what we 
will use to parametrize maximal states. In Ref. p4| , it is 
shown that, given a fixed set of eigenvalues, all states that 
maximize one of the three entanglement measures (en- 
tanglement of formation, negativity or relative entropy) 
automatically maximize the other two. It was further 
shown that the global unitary transformation that takes 
arbitrary states into maximal ones has the form 



U={U X ® U 2 )TD^ 



(4.1) 



where U\ and XJi are arbitary local unitary transforma- 
tions, 



/ 1\ 
^ 7§ °, 

V 1 0/ 



(4.2) 



is a unitary diagonal matrix, and $ is the uni- 
tary matrix that diagonalizes the density matrix p, i.e., 
p = $A<f>^, where A is a diagonal matrix, the diagonal 
elements of which are the four eigenvalues of p listed in 
nonincreasing order (Ai > A2 > A3 > A4). Hence, the 
general form of a density matrix that is maximal, given a 
set of eigenvalues, is (up to local unitary transformations) 




/A 4 

Q Ai+A 




2 

A1-A3 





Al-A; 

2 

Al+A; 

2 




\ 




a 2 ; 



(4.3) 



This matrix is locally equivalent to the form 




Q 



§ 

a 

6 

£ + 



(4.4) 



with x+j = (Ai+Aa)/2, r — X\— A3, a — A2, and b — A4. The 



above derivation justifies the Ansatz form (4.5) used in 
Ref. | fl2| to derive the entanglement of formation vs. lin- 
ear entropy MEMS. We remark that one may as well use 
the four eigenvalues A^'s as the par ame trization. Nev- 
ertheless, the form (4.4), as well as (4.5), can be nicely 
viewed as a mixture of a Bell state \(j> + ) with some diag- 
onal separable mixed state. 
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B. Entanglement-versus-linear-entropy frontiers 

We begin by measuring mixedness in terms of the lin- 
ear entropy, and comparing the frontier states for various 
measures of entanglement. 



1. Entanglement of formation 

The characterization of physical states in terms of their 
entanglement of formation and linear entropy was intro- 
duced by Munro et al. in Ref. O]. (Strictly speaking, 
they considered the tangle rather than the equivalent en- 
tanglement of formation.) Here, we shall consider yet an- 
other equivalent quantity: concurrence (see Sec. II A). In 
order to find the frontier, Munro et al. proposed Ansatz 
states of the form 



/x + 



PAn 







V 



a 
6 
y + 



(4.5) 



where x, y,a,b,r > and x + y + a + b + r = 1. They 
found that, of these, the subset 



PmEMS:E f , 5 L — 



Pi(r) = 



1-r-O 





x 2 




forf 



< r < 1 
for < r < | 



Ph(t-) = 



(\ 











y 



(4.6a) 



(4.6b) 



lies on the boundary in the tangle vs. linear-entropy plane 
and, accordingly, named these MEMS, in the sense that 
these states have maximal tangle for a given linear en- 
tropy. We remark that at the crossing point of the two 
branches, r = 2/3, the density matrices on either side 
coincide. 

In Fig. [| we plot the entanglement of forma- 
tion/concurrence vs. linear entropy for the family of 
MEMS ( |4.6|) ; this gives the frontier curve. For the sake 
of comparison, we also give the curve associated with the 
family of Werner states of the form 




0.2 0.4 0.6 0.8 

Linear entropy 




0.2 0.4 0.6 0.8 1 

Linear entropy 

FIG. 3: Entanglement frontier. Upper panel: entanglement 
of formation vs. linear entropy. Lower panel: concurrence 
vs. linear entropy. The states on the boundary (solid curve) 
are PmemS;B f ,s l • A dot indicates a transition from one branch 
of MEMS to another. The dashed curve below the boundary 
contains Werner states. 



2. Relative entropy as the entanglement measure 



To find the frontier states for the relative entropy of 
entanglement we aga in turn our attention to the maximal 
density matrix (4.4). For this form of density matrix 
the linear entropy is given (with x expressed in terms of 
a,b,r) by 



S L = -(-3a 2 + 2a(l-&) + (l-6)(l + 36)-r 2 ). (4. 



Pw 



-II = 



(±£ 

o ^ o 

1-r 

4 

i±£ 



1-r 
4 


S 






r+r / 
4 / 



(4.7) 



Evidently, for a given value of linear entropy these MEMS 
(which we shall denote by {MEMS : Ep, Sl}) achieve the 
highest concurrence. As the tangle r and entanglement 
of formation E-p are monotonic functions of the concur- 
rence, Eq. (4.6) also gives the boundary curve for t hese 
measures. This raises an interesting question: Is (4.6) 
optimal for other measures of entanglement? 



To calculate the relative entropy of entanglement, we 
need to determine the closest separable state to ([O]). 
It is simpler to do this analysis via several cases. We be- 
gin by considering the Rank-2 and Rank-3 cases of (fl.4|). 
We set b = (A 4 = 0) and express x in terms of a and r 
in the density matrix, obtaining 



P = 




a 











I — a 
2 / 



(4.9) 
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and thus find that the closest separable density matrix 
a* is given by Jul 



C = 
D = 
E = 



<C D\ 

E 

1-2C-E 

vD C) 

(l+a) (l-a 2 -r 2 ) 
2(l+a-r)(l+a+r)' 

a (l + a) r 
(1+a-r) (1+a+r)' 

a (l+a) 2 
(1 + a-r) (1 + a+r)' 



(4.10a) 

(4.10b) 
(4.10c) 
(4.10d) 



The relative entropy of entanglement is now simply given 
by 



E R (p) = Tr (p log p - p log <t* 
l + a (l + a) 2 -r 2 7 



log 



2 ~° (l + a) 2 
with the linear entropy being given by 

5 L = ?(l + 2a-3a 2 -r 2 ) 



r 1+a+r 
- log 

2 5 1+a-r 



,(4.11) 



(4.12) 



subject to the constraint (a + r) < 1. For the Rank- 2 
case a = 1 — r (b — x — 0), and the r esult ing solution 
is the Rank-2 matrix pi(r) given in Eq. (4.6) with 1/2 < 
r < 1. We remark that this Rank-2 solution is always 
a candidate MEMS for the three entanglement measures 
that we consider in this Paper. In order to determine 
whether or in what range the Rank-2 solution achieves 
the global maximum, we need to compare it with the 
Rank-3 and Rank-4 solutions. 

By maximizing E R (p) for a given value of Sl, we find 
the following stationary condition: 



r log 



(l+a) 2 



(1 



(3a - 1) loj 



1 



1 



(4.13) 



Given a value of Sl, we can solve Eqs. (4.12) and (4.13), 
at least numerically, to obtain the parameters a and r, 
and hence, from Eq. (19|), the Rank-3 MEMS. However, 
if the constraint inequality a + r < 1 turns out to be 
violated, the solution is invalid. 

We now turn to the Rank-4 case. It is straightforward, 
if tedious, to show that the Werner states, Eq. (4.7), 
obey the stationarity conditions appropriate for Rank 4. 
However, it turns out that this solution is not maximal. 

To summarize, the frontier states, which we d enote by 
{MEMS : E R ,S L }, are states of the form (Q; the de- 
pendence of the parameters a and r on Sl is shown in 
Fig. H In Fig. ||, we show the resulting frontier, as well as 
curves corresponding to non-maximal stationary states. 
The frontier states have the following structure: (i) for 
S L < 0.5054 they are the Rank-2 MEMS of Eq. (O) but 
with r restricted to the range from 1 (at Sl = Ojto ap- 
proximately 0.7459 (at S L ~ 0.5054); (ii) for S L > 0.5054 



0.6 




0.4 0.6 0.8 1 

Linear entropy 

FIG. 4: Dependence of a and r of the frontier states on linear 
entropy. The dotted line indicates the transition between two 
branches of MEMS. 




0.4 0.6 

Linear entropy 

FIG. 5: Entanglement frontier: relative entropy of entan- 
glement vs. linear entropy. The frontier states (solid curve) 
are PmemS;E r ,s l - The dot indicates the transition between 
branches of MEMS. 



the MEMS are R ank 3, with parameters a and r satis- 
fying Eqs. (4.13) and ( 4.1 2| ) at each value of Sl, and 
(a, r) ranging between approximately (0.3056, 0.7459) (at 
S L ~ 0.5054) and (1/3,0) (at S L = 8/9). As noted pre- 
viously, beyond Sl = 8/9 there are no entangled states. 
As the inset of Fig. [l] shows, the parameter a can be re- 
garded as a continous function of parameter r £ [0, 1]. 
The two branches of the solution, (i) and (ii), cross at 
(Sl,E R ) ~ (0.5054,0.3422); at this point, the states on 
the two branches coincide, 



'0.372947 0.372947\ 

0.254106 



v0. 372947 0.372947/ 



(4.14) 



Just as in the case of entanglement of formation vs. linear 
entropy, the density matrix is continuous at the transition 
between branches. 

We remark that the curve generated by the states 
{MEMS: .Bp, S L }, when plotted on the E R vs. S L plane, 
falls just slightly below that generated by the states 
{MEMS : £ R ,S L } for S L > 0.5054 (and coincides for 
smaller values of Sl). We also remark that the parame- 
ter r turns out to be the concurrence C of the states, so 
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0.2 0.4 0.6 0. 

Linear entropy 

FIG. 6: Differe nce in entanglement (Ep — E-n) vs. Sl for the 
MEMS in Eq. (4.6) and Werner states. The solid curve shows 
states from pmemS;B f ,s l 1 the dashed curve shows the Werner 
states. 



that Fig. |I] can be interpreted a plot of the concurrence 
of the frontier states vs. their linear entropy. By compar- 
ing this concurrence vs. linear entropy curve to that in 
Fig. ||, we find that the former lies just slightly below the 
latter for Sl > 0.5054 (and the two coincide for smaller 
values of Sl), the maximal difference between the two 
being less than 10~ 2 . 

It is evident that, for a given linear entropy, the rela- 
tive entropies of entanglement for both {MEMS : Er, Sl} 
and {MEMS :Ep, Sl} are significantly less than the cor- 
responding entanglements of formation. In fact, for small 
degrees of impurity, the entanglements of formation for 
the two MEMS states are quite flat; however, the rela- 
tive entropies of entanglement fall quite rapidly. More 
specifically, for a change in linear entropy of ASl =0.1 
near Sl = we have that AEp w 0.05 (see Fig. |3|) and 
AEr w 0.2 (see Fig. ||). As the curves of the states 
{MEMS : E F , S L } and {MEMS : E R , S L } are very close 
on the two planes, E-p vs. Sl and E R vs. Sl, we show, 
in Fig. H the entanglement difference Ep — En for the 
states {MEMS : £"f,Sl}, and compare it with the cor- 
responding difference for the Werner states. While it is 
clear that Epi(p) < Ep(p), for certain values of the linear 
entropy the difference turns out to be quite large, this 
difference being uniformly larger for {MEMS : Ep,Si,} 
than for the Werner state; see Fig. This raises sev- 
eral interesting questions, e.g., whether, for a given value 
of Ep, the undistillable entanglement (the difference be- 
tween the entanglement of formation and the relative en- 
tropy of entanglement is the lower bound) increases as 
states become more mixed. 

As we have seen, Werner states are not frontier states 
either in the case of entanglement of formation or in the 
case of relative entropy of entanglement. By contrast, 
as we shall see in the next section, if we measure en- 
tanglement via negativity, then for a given amount of 
linear entropy, the Werner states (as well as another 
Rank-3 class of states) achieve the largest value of en- 
tanglement. Said equivalently, the Werner states belong 
to {MEMS:N,S L }. 



3. Negativity 

In order to derive the form of the MEMS in the case 
of negativity, we again consider the density matrix of the 
form (4.4), for which it is straightforward to show that 



the negativity A is given by 



A = max{0, y/(a - b) 2 + r 2 - (a + b)}. (4.15) 

Furthermore, because we aim to find the entanglement 
frontier, we can simply restrict our attention to states 
satisfying A > 0, i.e., to states that are entangled pi] . 
Then, by making A stationary at fixed Sl and with the 
constraint 2x + a + b + r= 1, we find two one-parameter 
families of stationary states (in addition to the Rank- 
2 MEMS, which are common to all three entanglement 
measures). The parameters of the first family obey 



1 - 4x. 



(4.16) 



When expressed in terms of parameter r, the density ma- 
trix takes the form 



/AlEMS:iV,S L 



/ 1+r 
' 4 






V 





1-r 

4 






o § \ 


o ifV 



(4.17) 



which are precisely the Werner states in Eq. (4.7). For 
the second solution, the parameters obey 



4-2V3^2+T , n l + VS^+l r 
a= ,b = 0,x = - . (4.18) 



G 



G 



When expressed in terms of parameter r, the density ma- 
trix takes the form 



I 6 







(2) 

PmEMS:N,St 



V 



4-2^3^+1 
6 






2 




1 + ^3^+1 



(4.19) 



We remark that the two solutions give the same bound on 
the negativity for a given value of linear entropy. The re- 
sulting frontier in the negativity vs. linear entropy plane 
is shown in Fig. |?]. We have shown by numerical ex- 
ploration [^7j that the family states found in the previ- 
ous paragraph do indeed achieve the maximal negativity. 
The results of this exploration are indicated by points in 
the inset in Fig. ^, and compelling evidence for the max- 
imally of the negativity is furnished by the fact that no 
points lie beyond the frontier. 

Thus, the states {MEMS: A, S L } on the boundary in- 
clude, up to lo cal u nitary transformatio ns, bo th Werner 
states in Eq. ( |4.17| ) and states in Eq. (|4.19| ). We also 
plot in Fig. g the curve belonging to {MEMS : Ep, S L }; 
note that it falls slightly below the curve associated with 
{MEMS : A, Sl} and that it has a cusp, due to the struc- 
ture of t he states, at the value 2/3 for the parameter r in 
Eq. ( [4.61 ). Here, we see that maximally entangled mixed 
states change their form when we adopt a different en- 
tanglement measure. 
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Linear entropy 

FIG. 7: Entanglement frontier: negativity vs. linear en- 
tropy. States on the boundary (full line) are Pmems at s l anc ' 

PMEMS:iv,s L - Tne dashed line comprises PmemS:E f ,s l • In- 
set: the randomly generated states confirm the location of 
the frontier. 



Entanglement- versus- von-Neumann-entropy 
frontiers 



We continue this section by choosing to measure 
mixedness in terms of the von Neumann entropy, and 
comparing the frontier states for various measures of en- 
tanglement. 



1. Entanglement of formation 

To find this frontier, we consider states of the 
form (4.4), and compute for them the concurrence and 
the von Neumann entropy: 



C=r-2Vab, (4.20a) 
Sy = —a log a — b \ogb—x log x — (x+r) log(x+r). (4.20b) 

Note that the parameters obey the normalization con- 
straint 2x+a + b+r = 1. 

As we remarked previously, the Rank 2 MEMS is al- 
ways a candi date. For the Rank-3 case, we can set 6 = 
in Eq. ( 4.20 ). By maximizing C at fixed Sy, we find a 
stationary solution: 

(i) r = C, x = (4-3C-v / 4-3C* 2 )/6, and a = (V4-3C 2 - 
l)/3; the resulting density matrix is 



Pi = 



2 



1 

VEME-l o 


4-V4-3C 2 
6 



V T 



o 



(4.21) 



For the Rank-4 case (b ^ 0), the stationarity condition 
can be shown to be 



C 

O 0.8 
'«*— < 
(0 

E 0.6 



h- 0.4 

o 



C 0.2 
LU 



Rank3 



Rank2 
Transcendental 



Werner 



0.2 0.4 0.6 0.8 

von Neumann entropy 




0.2 0.4 0.6 0.8 1 

von Neumann Entropy 

FIG. 8: Entanglement frontiers. Upper panel: entanglement 
of formation vs. von Neumann entropy. Lower panel: con- 
currence vs. von Neumann entropy. The branch structure is 
described in the text. 



where u 



\f a/ (x+r), v — y x/(x+r), and w = 
yjb/ (x + r). There are two solutions, due to the two- 
to-one property of the function z\ogz for z £ (0,1). The 
first one is (u = v = w). 

(ii) a = b = x = (l-C)/6, and r = (l+2C)/3, which can 
readily be seen to be a Werner state as in Eq. (4.7) or, 
equivalently, 



pa 



I 2+C 
' 6 




1+2C 





l-C 
6 









l-C 
6 





1+2C \ 
6 » 




2+C 
6 



(4.23) 



Being the concurrence, C is restricted to the interval 
[0,1]. The second solution is transcendental, but can 
be solved numerically. 

In Fig. H we compare the four possible candidate solu- 
tions, and find that the global maximum is composed of 
only (i) and (ii). We summarize the states at the frontier 
as follows: 



ulog(u) = u;log(w), (4.22a) 
2u\og(u) = (u + w)log(v), (4.22b) 



PMEMS:E f ,Sv 



for < C < C*; 
for C* < C < 1. 



(4.24) 
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Note the crossing point at (C,Sy) = (C*,Sy(C*)), at 
which extremality is exchanged, so the true frontier con- 
sists of two branches. It is readily seen that C* is the 
solution of the equation Sy(p\(C)) = Sy(pn(C)), and 
the approximate numerical values of C* and the corre- 
sponding S v are 0.305 and 0.741, respectively. 

The resulting form of MEMS states is peculiar, in 
that, even at the crossing point of two branches on the 
entanglement-mixedness plane, the forms of matrices on 
the two branches are not equivalent (one is Rank 3, the 
other Rank 4). This is in contrast to the {MEMS : 
Ef,Sl}. This peculiarity can be partially understood 
from the plot of the two mixedness measures, Fig. || as 
the value of the von Neumann entropy rises, there are 
fewer and fewer Rank-3 entangled states, and above some 
threshold, no more Rank-3 states exist, let alone entan- 
gled Rank-3 states. There are, however, still entangled 
states of Rank 4. Hence, if Rank-3 states attain higher 
entanglement than Rank-4 states do when the entropy 
is low, a transition must occur between MEMS states of 
Rank 3 and Rank 4. 

From Fig. || it is evident that beyond a certain value of 
the von Neumann entropy no entangled states exist. This 
value can b e readily obtained by considering the MEMS 
state Q4.23D at C = 0, 







? *\ 

| 

o o I 



(4.25) 



for which Sy = -(1/2) log 4 (l/12) » 0.896. 

As an aside, we mention a tantalizing but not yet fully 
developed analogy with thermodynamics pq] . In this 
analogy, one associates entanglement with energy and 
von Neumann entropy with entropy, and it is therefore 
tempting to regard the MEMS just derived as the analog 
of thermodynamic equilibrium states. If we apply the 
Jaynes principle to an ensemble in equilibrium with a 
given amount of entanglement then the most probable 
states are those MEMS shown above. 




FIG. 9: Dependence of _Er, Sy and a on r for the Rank-3 
maximal states. 
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von Neumann entropy 

FIG. 10: Entanglement frontier: relative entropy of entan- 
glement vs. von Neumann entropy. The solid curve is the 
frontier. The branch structure is described in the text. 



rescaling of the constraint-enforcing Lagrange multiplier. 
Thus, in maximizing at fixed Sy, we arrive at the 
stationarity condition 



log 



(l+a) : 



log 



1— a— r 
l—a+r 



, 1+a+r (l-a) 2 -r 2 
dog— log- 



1+a—r 



4a 2 



(4.27) 



2. Relative entropy of entanglement 

Let us now find the frontier states for the case of rela- 
tive entropy of entangleme nt . To do this, we first consider 
the Rank-3 states in E q. (|4.9| ), for which the relative en- 
tropy is given by Eq. ( 4. lip . For these states, the von 
Neumann entropy is given by 



l—a+r l—a+r 
s v = o — lo S — s — 



a log a 



1— a— r 1— a— r 



(4.26) 



where the log function is taken to be base four. Even 
though the log functions in Sy and use different bases, 
the stationary condition for the parameters r and a does 
not change, because the difference can be absorbed by a 



We can solve for the parameter a as a function of r e 
[0, 1], at least numerically; the result is shown in Fig. ||, 
along with Sy and 1?r. 

Turning to the Rank-4 case, it is straightforward, if 
tedious, to show that the Werner states satisfy the cor- 
responding stationarity conditions. In order to ascertain 
which rank gives the MEMS for a given Sy, we compare 
the stationary states of Rank 2, 3, and 4 in Fig. [l0| Thus, 
we see that for Sy < Sy ~ 0.672 the frontier states are 
given by the Rank-3 states, whereas for Sy > Sy the 
frontier states are given by the Werner states ([O]) with 
the parameter r ranging from approximately 0.6059 down 
to 0. At the crossing point, (Sy, E&) ~ (0.672, 0.124), the 
MEMS undergo a discontinuous transition; recall that we 
encountered a similar phenomenon in the case of entan- 
glement of formation vs. von Neumann entropy. 
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FIG. 11: Entanglement frontier: negativity vs. von Neumann 
entropy. The solid curve is the frontier. The broken curve 
represents the Rank-2 candidate states. 



3. Negativity 



We saw in Sec. [VB 3 that there is a pair of families of 
MEMS which differ in rank but give the identical frontier 
in the N vs. Sx plane. It is interesting to see what hap- 
pens for the combination of negativity and von Neumann 
entropy. 

Once again, we begin with states of form (4.4), for 



which the neg ativi ty an d the v on Neumann entropy are 
given in Eqs. ( |4.15| ) and ( 4.20b ), respectively. By making 
N stationary at fixed Sv, we are able to find only one 
solution (in addition to the Rank-2 candidate): a = b = 
x. Expressing the resulting density matrix, as we may, in 
terms of the single parameter r, we arrive at the following 
candidate for the frontier states: 



PMEMS:N,S\ 



/ l+r 
' 4 








1-r 
4 









1-r 
4 





§\ 




l+r 1 
4 / 



(4.28) 



where < r < 1, i.e., the Werner states. 

The resulting frontier in the negativity vs. von- 
Neumann-entropy plane is shown in Fig. [n] which, for 
comparison, also shows the curve for the Rank-2 candi- 
date. 



V. CONCLUDING REMARKS 

In this Paper we have determined families of maximally 
entangled mixed states (MEMS, i.e., frontier states, 
which possess the maximum amount of entanglement for 
a given degree of mixedness). These states may be use- 
ful in quantum information processing in the presence 
of noise, as they have the maximum amount of entan- 
glement possible for a given mixedness. We considered 
various measures of entanglement (entanglement of for- 
mation, relative entropy, and negativity) and mixedness 
(linear entropy and von Neumann entropy). 

We found that the form of the MEMS depends heavily 
on the measures used. Certain classes of frontier states 




0.2 0.4 0.6 0.8 

Linear entropy 
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von Neumann entropy 

FIG. 12: Violation of Bell's inequality for various fami- 
lies of states: (a) {MEMS : E f ,Sl}, (b) Werner states, 



( c ) Pmems.n.s > (d) pi i n E q. (4.21), and (e) the Rank-2 Bell 



diagonal states in Eq. (5.3) 



(such as those arising with either entanglement of forma- 
tion or relative entropy of entanglement vs. the von Neu- 
mann entropy) behave discontinously at a specific point 
on the entanglement-mixedness frontier. Under most of 
the settings considered, we have been able to explicitly 
derive analytical forms for the frontier states. 

For entanglement of formation and relative entropy- 
and for most values of mixedness, we have found that the 
Rank-2 and Rank-3 MEMS have more entanglement than 
Werner states do. On the other hand, at fixed entropy no 
states have higher negativity than Werner states do. At 
small amounts of mixedness, the {MEMS:^f, <Sl} states 
"lose" entanglement with increasing mixedness at a sub- 
stantially lower rate than do the Werner states. However, 
when the entanglement is measured by the relative en- 
tropy, the difference in loss-rate is significantly smaller. 

From Eq. ( [2.1C ) it is tempting to assert that in the 
case of LS entanglement vs. mixedness the frontier states 
should be the same as those in the case of entanglement 
of formation vs. mixedness. However, as we do not know 



whether Eq. (4.5) (up to local unitary transformations) 
exhausts all maximal states for LS entanglement, further 
investigation of this point is needed. 

Having characterized the MEMS for various measures, 
it is worthwhile considering them from the perspective 
of Bell-inequality violations. To quantify the violation of 
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Bell's inequality, it is useful to consider the quantity 
B = max {E(a, b) + E(a, V) + E(a! , b) - E(a' , b') } , 

a, a', 6,6' 

(5.1) 

where E(a, b) = (a ■ a® a -b) and the vectors a and a! (b 
and 6') are two different measuring apparatus settings for 
observer A (observer B). If B > 2 then the corresponding 
state viol ates Bell's inequality. For the density matrix of 
the form ( |4.4| ) it is straightforward J29| to show that the 
quantity B is given by 



B 



(5.2) 



In Fig. [12], we plot B vs. linear and von Neumann en- 
tropies for several families of frontier states. As a com- 
parison, we also draw the violation by the following 
Rank-2 state (which is diagonal in the Bell basis): 

PB = r\cP+)(c(>+\ + (l-r)\<p-){<t>-\, re [0,1]. (5.3) 

This state, although not belonging to any of the fam- 
ilies of frontier states derived previously, turns out to 
achieve the maximum possible violation, as a function 
of linear entropy. On the other hand, the Werner states 
appear to achieve maximal violation in the case of von 
Neumann entropy. Ekert's application of Bell's inequali- 
ties to quantum cryptography [|| , together with the dis- 
cussions of the present paragraph, suggests that MEMS 
may be relevant to quantum communication. 

Another natural application for which entanglement is 
known to be a critical resource is quantum teleportation. 
How do these frontier MEMS teleport, compared with the 
Werner and Rank-2 Bell diagonal states? If we restrict 
our attention to high purity situations (i.e., to states with 
only a small amount of mixedness) then it is straight- 
forward to show that, e.g., {MEMS : E-p, S^} states tele- 
port average states better than the Werner states do, but 
worse than the Rank-2 Bell diagonal state does. Part of 
the explanation for this behavior is that standard tele- 
portation is optimized for using Bell states as its core 
resource. 

It is also interesting to note that for certain combina- 
tions of entanglement and mixedness measures, as well 
as the Bell inequality violation, the Rank-2 candidates 
fail to furnish MEMS. Thus, these states seem to be 
less useful than other MEMS. However, from the per- 
spective of distillation, these states are exactly quasi- 
distillable |30[ and can be useful in the presence of 
noise because they can be easily distilled into Bell states. 
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APPENDIX A: NUMBER OF NEGATIVE 
EIGENVALUES OF THE PARTIAL TRANSPOSE 

OF p 

In this appendix we address the result that for C 2 £g> C 2 
systems the partial transpose of any density matrix p has 
at most one negative eigenvalue. In fact, we shall consider 
the result from two perspectives. 

First, we build upon results (Theorem 3, in particular) 
contained in Ref . |}l| , from which it follows that it is suf- 
ficient to consider (i) Bell diagonal states and (ii) states 
of the form 



/a + c d 


b-c 

V d a-b) 



(Al) 



For the latter case, straightforward calculation shows 
that the partially transposed matrix can have one neg- 
ative eigenvalue when d =/= and that it does not 
have negative eigenvalues when d — 0. For the for- 
mer case, suppose that the Bell diagonal state has the 
four eigenvalues Ai, A2, A3, A4, in nonincreasing or- 
der. Then it is straightforward to see that the corre- 
sponding partially transposed matrix has the four eigen- 
values (Ai + A 2 + A 3 - A 4 )/2, (Ai + A 2 - A 3 + A 4 )/2, 
(Ai - A 2 + A 3 + A 4 )/2, (-Ai + A 2 + A 3 + A 4 )/2. Thus, it 
can have at most one negative eigenvalue. 

A second perspective is provided by first invoking the 
LS decomposition ( |2.8| ) and then making a Schmidt de- 
composition, via the local unitary transformation U a <8> 
Ub, of the pure entangled part \ip e )- In this way, the pure 
part becomes 



c/A^c/slV'e) = Vp|oo) + v / i 3 p|ii), (A2) 

where 1/2 < p < 1. Meanwhile, p is transformed into 



p' = (U A ® U B ) P [U\ ® U%) 

( V O0vpFri\ 





Ap' s + (1-A) 



, (A3) 



Vp(i-p)oo 1-p ) 

S U\®U\ is still separat 
its partial transpose remains positive semi-definite 



where p' s = Ua ® Ub p s U\ ® W B is still separable, and 
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Taking, then, the partial transpose, the density matrix 
is transformed into 



P = *Ps 



(1-A) 



(V ( 
C 

\o o 










_° \ 

-p) o 



l-pj 



,(A4) 



and we no te that the last matrix has ei genvalue s Ai = 
Pi A 2 = VpQ^-p): A 3 = A4 = in (as 

P > 1/2) nonincreasing order. As the Hermitian matrix 
p' s TB retains positive-semi-definiteness, we can employ a 
well-known result in matrix analysis |32|] that for n x n 



Hermitian matrices A and B, with i? being positive semi- 
definite, the eigenvalues of (A+B) and A, when arranged 
in non-ascending order, obey 



X k (A + B)>X k {A) fork = 1,2,. 



(A5) 



Hence, identifying B with Xp' s B and A with the product 
of (1— A) and the matrix in Eq. (A4), we immediately see 



that Xi(j/ Ta ) > A 2 ( P ' Tb ) > X 3 ((F B ) > (l-A)(l-p) > 
and, thus, p' Tb (or equivalently p Tf3 ) can have at most 
one negative eigenvalue. Thus, the negativity for C 2 <g)C 2 
systems can then be written as N = 2 max{0, — A^p 1 ^ )}. 
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